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Abstract 


As one class of the most important objects in the universe, magnetars can produce a lot of different frequency 
bursts including X-ray bursts. In Cai et al., 75 X-ray bursts produced by magnetar SGR J1935+2154 during an 
active period in 2020 are published, including the duration and net photon counts of each burst, and waiting time 
based on the trigger time difference. In this paper, we utilize the power-law model, dN (x)/dx c (x + xo) ^5, to fit 
the cumulative distributions of these parameters. It can be found that all the cumulative distributions can be well 
fitted, which can be interpreted by a self-organizing criticality theory. Furthermore, we check whether this 
phenomenon still exists in different energy bands and find that there is no obvious evolution. These findings further 
confirm that the X-ray bursts from magnetars are likely to be generated by some self-organizing critical process, 
which can be explained by a possible magnetic reconnection scenario in magnetars. 
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1. Introduction 


Magnetars are a group of isolated neutron stars, which have 
extremely powerful magnetic fields (Thompson & Dun- 
can 1995). The dissipation of the magnetic field provides 
energy for magnetars. They are characterized by strong 
variability on several timescales and exhibit large variations 
across the electromagnetic spectrum, especially in X-ray and 
soft gamma-ray energy bands, ranging from a few milliseconds 
to a month (Kaspi & Beloborodov 2017). SGR J1935+2154 
was discovered in 2014 when Swift-BAT (Burst Alert 
Telescope) was triggered by a short burst from Galactic plane 
(Stamatikos et al. 2014). Based on a full multi-wavelength 
radio study of continuous and persistent emission of SGR 
J1935+2154 reported by Kothes et al. (2018), Gaensler (2014) 
found that it was associated with supernova remnant SNR 
G57.24-0.8. Zhou et al. (2020) corrected the distance from 
SGR J1935+2154 to SNR G57.2+0.8 and found that the 
distance between them may only be 6.6+0.7 kpc, which is 
closer than the previous hypothesis of 10 kpc. 

Since 2020 April 27, SGR J1935+2154 entered a new 
period of activity, in which the magnetar produced several 
outbursts, including fast radio burst (FRB) 200428 (Bochenek 
et al. 2020; CHIME/FRB Collaboration et al. 2020). The 
observation of FRB 200428 first confirmed the prediction that 
magnetar may be one of the candidates of the origin of FRBs 


(Bailes et al. 2021; Cooper et al. 2022). In the days and months 
following FRB 200 428, several telescopes made sustained 
observations. However, there is no significant single radio 
pulse (Bailes et al. 2021). Notably, the Insight-HXMT satellite 
detected a non-thermal X-ray burst associated with FRB 
200 428 and then it was identified as the emission from SGR 
J19354-2154 (Li et al. 2021). Moreover, Insight-HXMT has 
detected 75 X-ray bursts (1-250 keV) from SGR J1935+2154 
after FRB 200428 and found that one of them shows a similar 
peak spectral energy from the X-ray burst associated with FRB 
200428 (Cai et al. 2022a, 2022b). They found that the 
cumulative distribution of the fluence of all 75 bursts can be 
well fitted by a power-law with an index of 0.764 + 0.004, 
which is consistent with the result of SGR J19354-2154 bursts 
by Fermi/GBM reported by Lin et al. (2020). 

Studies show that X-ray events often occur in some high-energy 
celestial burst events, e.g., X-ray flares from X-ray binaries or other 
systems, especially in the nearest star, the Sun (Aschwanden 2011; 
Shibata & Magara 2011), the most violent celestial explosions 
known in the universe: gamma-ray bursts (GRBs; Burrows et al. 
2005; Falcone et al. 2006; Zhang et al. 2006; Mu et al. 2016; Yi 
et al. 2016, 2017; Tang et al. 2021; Li et al. 2023), some active 
galactic nuclei events (AGNs; Rees 1984; Yan et al. 2018), and the 
tidal disruption events (TDE), such as Swift J1644+57 (Bloom 
et al. 2011; Burrows et al. 2011; Zheng et al. 2020). 
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As is well known that the solar X-ray flares could be originated 
from the magnetic reconnection process, Wang & Dai (2013) 
found that GRB X-ray flares show the similar power-law 
distributions like solar flares in terms of the waiting time, energy 
and duration, respectively. It is believed that X-ray flares from 
GRBs and solar flares may also have the same physical 
mechanism. Because they can also be estimated in a self-organized 
criticality (SOC) system (Bak et al. 1987; Bak & Tang 1989; Lu & 
Hamilton 1991; Aschwanden 2011; Shibata & Magara 2011). 

In addition, the same power-law distributions of X-ray flares 
or bursts have been found in other systems, for example, some 
black hole binary systems (Wang et al. 2015; Yan et al. 2018), 
type I X-ray bursts originated from X-ray binary systems with 
low-mass (Wang et al. 2017), some repeating FRBs (Wang 
et al. 2017; Zhang et al. 2020; Wang et al. 2023), and soft 
gamma repeaters (e.g., Cheng et al. 1996; Göğüs et al. 
1999, 2000; Cheng et al. 2020). The similar distributions for 
parameters of X-ray flares from supergiant fast X-ray transients 
are also reported by Sidoli et al. (2016) and Zhang et al. (2022). 
Their studies indicated that the power-law characters for the 
supergiant fast X-ray transients are important evidence of SOC. 

The SOC behaviors of high energy celestial bodies indicate 
that it is necessary to systematically analyze the properties to 
further study their physical origin. In this work, we study the 
statistical properties of 75 bursts generated from SGR J1935 
+2154 reported by Cai et al. (2022a), including the durations, 
waiting times, and net photon counts. In Section 2, we present 
the sample selection and data analysis methods. Our main 
results and discussion are shown in Section 3. Finally, our 
conclusions are given in Section 4. 


2. Data and Statistical Analyses 


From 2020-04-28T07:14:51 UTC to 2020-06-01T00:00:01 
UTC, Insight-HXMT conducted a 33 days ToO observation of 
SGR J1935-F-2154 (total span 2851.2 ks), during which several 
hundred short X-ray bursts were emitted from SGR J1935 
+2154. The detailed observation time is listed in Table 4 
reported by Cai et al. (2022a). They provide a comprehensive 
monitoring of the evolution of burst activity from SGR J1935 
+2154 with high temporal resolution and sensitivity over a 
very broad energy range (1—250 keV). In this work, we study 
three parameters, including duration, waiting time and net 
photon counts. The waiting time is defined as the difference 
between the beginning time of the i+ 1th and ith burst, that is, 
T waiting = PstartLi--l — Tar (Yi et al. 2016). 

In general, the differential distribution could be described by 
a threshold power-law distribution with the following Equation 


A GL prank eee (1) 
dx 


The cumulative distribution of such explosive events can be 
written as the integral of the whole number of events exceeding 
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a given value x, so the cumulative distribution function 
corresponding to Equation (1) can be expressed as (a, z 1) 
(Aschwanden 2015) 


Nam >x) = 1 + (Nas — D 


(x2 + xo)! ^: — ( + xo)! ^* 


(2) 


where o; is the variable parameter referring to the power-law 
index of the distribution, xo is a constant introduced by taking 
into account the threshold effect (e.g., incomplete sampling 
below xo, background contamination), Neny is the total number 
of events, x, and x» are the minimum and maximum values of 
x, respectively. The uncertainty of the cumulative distribution 
in a given bin i is approximately calculated as Ogum,i = 4 Mum. 
where Neum, is the number of events in the ith bin. 

The standard reduced chi-square (x2) goodness is used to 
confirm the best fit. The x, can be expressed as 


x As M2 
[Noum,th (xi) ; Neum,obs (xi)] (3) 


1 
Xv,cum m 
(ny ul Npar) i=l O cum,i 


for the cumulative distribution function (Aschwanden 2015), 
where n, is the number of logarithmic bins, n4, is the number 
of the free parameters, Neum,ops(x;) is the observed values, 
Neum,th(x;) is the corresponding theoretical values for cumula- 
tive distribution, respectively. It should be noted that the points 
below the threshold xo are simply noise and do not contribute to 
the accuracy of the best fitting of the power-law exponentials, 
so they are ignored when calculating the reduced chi-square. 

Owing to the limited number of bursts, we only analyze a 
cumulative distribution rather than a differential one. Only two 
free parameters appeared in the cumulative distribution 
function. In general, the cumulative distributions can be 
generated exceed the threshold xy due to incomplete sampling 
of the selected samples. Taking the threshold xo as the free 
parameter and adding the exponent a, to fit the cumulative 
distribution, the power-law exponent a, of the samples can be 
well constrained. At the same time, it should be indicated that 
the cumulative number distribution of Equation (1) is a power- 
law function with an exponent of a,, so it is an important quest 
of this work to deduce the power-law exponent from the data of 
these burst values. 

In this work, the python module pymce’ is utilized to compile 
the selected data and the Monte Carlo Markov chain (MCMC) 
method is used to obtain the confidence intervals of the fitting 
parameters. Due to oversimplified sampling at a low value 
threshold, the distributions of the selected parameters usually 
show a shallow part or a gap before the threshold xo 
(Aschwanden 2015; Wang et al. 2015). 


7 $ 2 
https: //pypi.org /project/pymc/ 
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3. Results and Discussion 


The scale-free power-law distribution of various events for 
different parameters, such as duration, energy, or luminosity, is 
one of the obvious characteristics of SOC systems (Paizis & 
Sidoli 2014; Zhang et al. 2022). According to the theoretical 
framework proposed by Aschwanden (2012), the concept of the 
fractal dimensions of an SOC avalanche system can be 
quantitatively connected to the cumulative frequency distribu- 
tions. In this framework, the power law index of SOC 
cumulative frequency distribution for different parameters is 
theoretically associated with the Euclidean space dimension 
S=1, 2, 3. 

Aschwanden (2012) reported theoretical indices such as the 
duration frequency distribution (œr) and the distribution of 
energy (ag), which aligning with ar = 2I and ag — I ; 
where S= 1, 2, and 3 are the Euclidean dimensions. We can 
see clearly that the indices are ay=ag=1 for S=1, 
ag — 1.29 and ar= 1.5 for S=2, ar=2 and ag=1.5 for 
S — 3, respectively. Here, the net counts reflect the amount of 
released energy of the X-ray bursts. 

In order to check the possible evolution across different 
energy bands, the distributions of three parameters in HE, ME, 
LE and the total energy bands are analyzed and shown in 
Figure 1. The detailed fitting results are listed in Table 1. As 
shown in Table 1 and Figure 1, the best-fitting power-law 
indexes of the cumulative distribution for the net photon counts 
in the total energy band, HE, ME and LE are 1.63 + 0.08, 
1.57 + 0.03, 1.58 + 0.03, and 1.58 + 0.03, which are basically 
consistent with the theoretical index ag = 1.5 for $ —3. The 
mean value of these indexes of four individual channels is 
1.59 + 0.02, which is larger than those of the fluence reported 
by Cai et al. (2022b). The steeper index of the net photon count 
distributions may be due to the difference of the range of the 
X-ray burst fluence. Similarly, the results of the duration and 
the waiting time in the four energy bands are 1.99 + 0.08, 
2.02+0.18, 1.99 0.11, 2.07 2:0.15, and 1.95+0.08, 
1.95 + 0.08, 1.95 + 0.07, 1.97 + 0.10, and also well consistent 
with the theoretical index og — 2 for S — 3. 

In addition, it can be found that there is no obvious power- 
law index evolution across different energy bands, which gives 
sufficient evidence to the SOC behavior in X-ray bursts 
produced by SGR J1935-F-2154. Moreover, the Euclidean space 
dimension S — 3 can be obtained according to the prediction of 
the FD-SOC model. It is the key point that the relationship 
among the indices of power-law distributions of different 
parameters of X-ray bursts is mainly dependent on the 
nonlinear scaling law among the self-organized critical 
parameters. 

Many previous studies have been performed on the statistical 
properties of SGRs, most of which focus on the distribution of 
burst energy (e.g., Cheng et al. 1996; Göğüs et al. 1999, 2000; 
Cheng et al. 2020). For example, Cheng et al. 1996 found that 
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SGR events and earthquakes share four unique statistical 
properties, including but not limited to power-law energy 
distribution and log-symmetric waiting time distribution. These 
statistical similarities suggest that SGRs should be powered by 
star-quakes like earthquakes. 

Göğüs et al. (1999) conducted a similar study of the bursts 
from SGR 1900+14 during its 1998-1999 active phase. They 
found that the distribution of fluence or energy of SGR 1900 
+14 aligns with a power law index of 1.66 and first reported 
the SOC behaviors in SGR bursts. Then, they further presented 
the statistics of bursts from SGR 1806-20 detected by the Rossi 
X-Ray Timing Explorer/Proportional Counter Array, the Burst 
and Transient Source Experiment, and the International 
Cometary Explorer (Göğüs et al. 2000). They found that the 
distribution of bursts' fluence observed with each instrument 
are well characterized by power laws with the indices 1.43, 
1.76, and 1.67, respectively. They proposed the hypothesis that 
the source of energy for the SGR bursts may not be any 
accretion or nuclear power, but rather crustquakes caused by 
the neutron star's evolving strong magnetic field. 

Similarly, the total counts or energy of SGR bursts are 
verified to have power-law-like size distributions. Cheng et al. 
(2020) found that the energy distribution of magnetar bursts 
can be well described by power-law functions with exponents 
of 1.84, 1.68, and 1.65 for the three events of SGR 71550-5418, 
SGR 1806-20, and SGR 1900+14, and the duration distribu- 
tions of them also show power-law forms with exponents of 
1.69, 1.72 and 1.82, respectively. Meanwhile, the distribution 
of waiting time can be described by a non-stationary Poisson 
process in which the occurrence rate increases exponentially. In 
Tables 2 and 3, we present a comparison of present relevant 
studies about SGRs, including the results of SGR J1935--2154 
and the other SGRs detected by different instruments. 

Lu et al. (2023) studied the energy spectra of these 75 bursts 
using a variety of models and found that the spectra of 
magnetar bursts are complex and diverse, and different types of 
bursts exhibit unique characteristics on the phase distributions 
due to different generation mechanisms. However, they found 
that most of the optimal energy spectrum can be fitted by the 
cutoff power-law and power-law models, which coincide with 
the power-law distribution on the statistical characteristics of 
the parameters. 

The SOC features have been found in FRB 121 102 (e.g., 
Cheng et al. 2020; Zhang et al. 2021). Recently, Wang et al. 
(2023) studied the frequency distributions of burst energy and 
waiting time of the two repeating FRBs 20121102A and 
20201124. According to the bimodal distributions of the waiting 
time, the bursts are divided into long and short parts. It is found 
that the two characteristics of both long and short bursts can be 
fitted by a power-law function and are well understood within 
the physical framework of an SOC system driven in a correlated 
way. They propose a possible trigger mechanism that favors the 
emission of the star's magnetosphere, the crustal failure of a 
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Figure 1. The cumulative distributions of the net count, duration and waiting time of the X-ray bursts from SGR J1935+2154. The gray region represents the 95% 
confidence level, the red line is the best fitting result, and the dashed line is marked as the threshold xo. Note that the numbers of data of each parameter adopted have 


been marked in each panel. 


The Best-fitting Power-law Indices (a) of These Parameters 


Table 1 
Band Net Counts 
Total(1-250 keV) 1.63 + 0.08 
HE(28-250 keV) 1.57 + 0.03 
ME(10-30 keV) 1.58 + 0.03 
LE(1-10 keV) 1.58 + 0.03 


Duration 
1.99 + 0.08 
2.02 + 0.18 
1.99+0.11 
2.07 + 0.15 


neutron star. The detection of FRB 200428 confirmed that at 
least some of the FRBs originate from magnetars (Bochenek 
et al. 2020; CHIME/FRB Collaboration et al. 2020). Wei et al. 
(2021) studied the properties of X-ray bursts from SGR J1935 
+2154 observed by the Gamma-ray Burst Monitor (GBM) on 
1.95 + 0.08 board and FRB 121 102 detected by the five-hundred-meter 
s ja M Aperture Spherical radio Telescope. They found that the high- 

energy components of FRB 121 102 and SGRs have a similar 


Waiting Time 
1.95 + 0.08 
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Table 2 
Results of Some Previous Studies of SGRs 
SGR Distribution? Function? ar QE References 
SGR 1627-41 cumulative PL 0.62 + 0.08 Woods et al. (1999) 
SGR J1550-5418 cumulative PL 1.69 + 0.03 1.84 + 0.03 Cheng et al. (2020) 
SGR 1806-20 differential PL m 1.6 Cheng et al. (1996) 
cumulative Threshold — PL 1.72 + 0.01 1.68 + 0.01 Cheng et al. (2020) 
SGR 1900+14 differential PL e 1.66 + 0.05 Göğüs et al. (1999) 
cumulative Threshold — PL 1.82 + 0.02 1.65 + 0.01 Cheng et al. (2020) 


Notes. 


* The purpose of distinguishing the forms of distributions is that because they are integral and differential to each other, the difference between the results of the same 


function fits is 1. 


> PLisa simple power law function (F « x~°*), and Threshold—PL is a threshold power law function of the two distribution forms used in this work. 


* The time-dependent parameters of a burst, including duration, waiting time, etc. 


The energy-related parameters of a burst, including energy, fluence, flux, photon-counts, etc. 


Table 3 
Comparison of Observation Results of the SGRs by Different Instruments 
SGR Distribution Function Og Instrument References 
SGR J19354-2154 differential PL 1.55 + 0.01 NICER Younes et al. (2020) 
cumulative PL 0.93 + 0.15 Fermi-GBM Rehan & Ibrahim (2023) 
cumulative Threshold—PL 1.63 + 0.08 Insight-HXMT this work 
SGR 1806-20 differential PL 1.43 + 0.06 PCA Göğüs et al. (2000) 
differential PL 1.76 + 0.17 BATSE 
differential PL 1.67 + 0.15 ICE 


scale-invariant behavior, which can be well interpreted by the 
same SOC framework with spatial dimension S = 3. 

Regarding the generation mechanism of magnetar X-ray 
bursts, Yuan et al. (2022) simulated the three-dimensional 
force-free electrodynamics of local Alfvén wave packets 
emitted by magnetar vibrations into the magnetosphere. They 
found that if the Alfvén wave packet propagates to a radius R 
and the total energy is greater than a certain value of 
magnetosphere energy, the wave would become very nonlinear 
and be ejected from the magnetosphere. The ejecta can carry a 
large energy of the initial Alfvén wave. Then it opens up the 
magnetospheric magnetic field lines, forming a current sheet 
behind them that connects back to the enclosed area. Magnetic 
reconnection will occur on these current sheets, which will 
cause plasma excitation and X-ray emission. We believe that 
the SOC process may occur during this period. 

In this work, we further studied the properties of X-ray 
bursts from SGR J1935+2154 observed by HXMT and found 
similar evidence of SOC power-law behaviors in all different 
energy channels. Therefore, we suggest that X-ray bursts from 
SGR J1935+2154 can be possibly considered as “avalanches” 
in SOC systems, in which the magnetic reconnection occurs 
during these X-ray bursts. 


4. Conclusions 


In this work, we have statistically analyzed three parameters 
of 75 X-ray bursts produced by SGR J1935+2154 during the 
active period beginning on 2020 April 27, including durations, 
waiting times, and net photon counts. We checked the 
cumulative distribution of these parameters in different energy 
bands. We found that all three parameters of X-ray bursts have 
similar power-law distributions, thus all can be explained by an 
SOC behavior. Moreover, it is found that there is no obvious 
power-law index evolution among different energy channels. In 
addition, the Euclidean space dimension S=3 has been 
obtained. Above all, we have obtained sufficient evidence that 
the X-ray bursts arise from a mechanism dominated by self- 
organizing critical systems, which will help us to further study 
the radiation mechanism of magnetars. 
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